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ABSTRACT 

Cosmological MHD simulations of galaxy cluster formation show a significant amplification of seed 
magnetic fields. We developed a novel method to decompose cluster magnetized turbulence into modes 
and showed that the fraction of the fast mode is fairly large, around 1/4 in terms of energy. This is 
larger than that was estimated before, which implies that cluster turbulence interacts with cosmic rays 
rather efficiently. We propose a framework to deal with electron and proton reacceleration in galaxy 
clusters that includes feedback on turbulence. In particular, we establish a new upper limit on proton 
and electron fluxes based on turbulence intensity. These findings, along with detailed modeling of 
reacceleration, will help to reconcile the observed giant radio haloes and the unobserved diffuse 7-ray 
emission from these clusters. 



1. INTRODUCTION 

Galaxy clusters are the largest virialized objects in the 
Universe. An important component of galaxy clusters 
is the intracluster meduim (ICM), a hot gas between 
galaxies which is fully ionized and contains magnetic 
fields and cosmic rays (CRs). A number of observa- 
tional techniques are used to detect various com ponents 
of th e ICM. Sunyaev-Zeldovich effect (e.g. ,Motl et al.l 
I2005f) and soft X-ray o bservations (e.g. see a review 
bv lMcNamara fc Nulsen..2007. ') are used to estimate hot 
gas temperature and density, radio observations of Fara- 
day rotation in the ICM measure magne t ic fields (e.g. 
Eiiek fc Owct|[2001 iBonafede et al.ll2010t iGovoni et al.l 
2006ll2010( ). Radio observations of Mpc scale diffuse ra- 
dio emissions, so called radio halos, indicate the pres- 
ence of cluster- wide rela t ivistic electrons in the ICM 
(e.g. ISch hckeiser et al.l f lQSTt ICarilU fc Tavloil [20021: 
Ferrari et al. 2008; Giovan niiiret al.ll2009t iFeretti et all 
2012tlvan Weeren et al.ll2012f )r 

The diffuse gamma-ray emission from the ICM, which 
should come from interactions of the proton compo- 
nent of CRs and the ga s, however, is not observed yet 
(|Ackermann et aLll2010[) . As the lower limits to diffuse 
7-ray emission are becoming better, it is harder to ex- 
plain both diffuse radio haloes and the lack of diffuse 
7-ray emission within a simple paradigm of the radio 
emission produced by secondary electrons (IReimer et ahl 
I2OOI iDoiI^ert et al.(l20l"lt iBrunetti fc Lazarianll2011bf ). 
Since the lifetimes of synchrotron electrons are much 
shorter than their mixing time in the cluster, the in- 
situ reacceleration mechanism is very desirable. One 
of such mechanisms being consi dered is a turbulent 
second -order acceleration (see, e.g. IBrunetti fc LazarianI 
201 lal). Radio relics s een in some clusters (see, e.g., 
van Weeren et all I2OIOI ) are thought to be associated 



with large-scale shocks. The explanation of radio haloes 
as shocks in projection, however, can be excl uded based 
on geometric arg uments ( Brunct ti et al.ll2008f) . or numer- 
ical simulations (iHoeft et al.i 120081: ISkillman et al.l [2011 
Vazza et al. 20il[). The second-order acceleration by tur- 
bulence, therefore, should be considered as a viable can- 
didate. The acceleration rate, which is proportional to 
(va/c)^ in the second order mechanism, was considered 
marginal in spiral galaxies, where 10 — 20km/s 

and the escape times from the relatively thin disk are 
short. In galaxy clusters, however, the Alfvenic speeds 
are higher va ^ lOOkm/s and the diffusion timescales are 
much longer due to their enormous size. 

The early theoretical model seeking to explain the 
spectr um of the diffuse radio halo of the Coma clus- 
ter by ISchhckeiser et al.l (|1987f) used the combination of 
first and second order Fermi acceleration and radiative 
synchrotron and inverse Compton losses to produce a 
volume-integrated frequency spectrum which could be 
approximated as oc ;y(3-r)/2 exp(— -y/ v/vc), with pa- 
rameters r and I'd defined in the aforementioned pa- 
per, depending on the details of the embedded shocks 
and the functional dependence of the diffusion coeffi- 
cients in energy. The above spectrum fits observations 
better than the single- and double-power laws ch arac- 
teristic of the primary (fja^fT977l iRephaelilTgr?) and 
secondary electron models (Jaffe 1977; Dennison 1980). 
The cutoff frequency depends on a number of parameters 
such as magnetic field, density, spatial diffusion coeffi- 
cient and the level of the background radiation field. Al- 
tho ugh a small fraction of galaxy clusters have radio ha- 
los (IFeretti et al.ll2012[). this fraction is higher for br ight 
clusters (iGio vannini et all ll999': Cass ano et al.|[2"007D . 

As a technique complementary to observations, full 
cosmological simulations of cluster formation and evolu- 
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Fig. 1. — A slice through the center of a simulated cluster data (3 Mpc across) - log density 



"icm^) (left), RMS magnetic field, 



/iG (middle) and RMS velocity, 10* cm/s (right). Typical sound speed Cs ~ 10* cm/s, typical Alfven speed va ~ lO'^ cm/s. 



tion, including mergers from infalling objects and accre- 
tion of gas, are becoming more and more popular recently 
(see, e.g. [N agai et al. 2007; Borgani & Kravtsov 2009; 
Dolag fc Stasvszvnll2009t IXu et all 120101 : iDonnert et all 



20111 , etc). These simulations produce magnetic fields 



which are roughly consistent with the observed rotation 
measures (see, e.g. iXu et al.|[20T2l ) . Thus, one can hope 
to elicit the properties of small-scale ICM turbulence and 
predict the effects of second-order acceleration. Obvi- 
ously, it is impossible to reach scales of CR gyroradius in 
a cosmological cluster simulation, so the subgrid model 
of turbulence has to be adopted. 

In §2 we discuss the origin of cluster magnetic field, in 
§3 we present the setup of our cosmological simulations, 
in §4 we describe compressible turbulence in the ICM, in 
§5 we propose a novel technique to decompose ICM tur- 
bulent perturbations into MHD modes, in §6 we discuss 
diffusive CR acceleration by the fast mode and in §7 we 
summarize our findings. 

2. SMALL-SCALE DYNAMO 

The strength of primordial magnetic fields is still 
unknowiQ with the upper bound around 10~ ^ G 
(jPaoletti fc Finelli|[20TTI; [Schleicher fc MiniatillMll) . As- 
suming that clusters formed as a result of a gravitational 
collapse of an unmagnetized gas, a proper mechanism 
of generating observed magnetic fields in the ICM has 
to be pointed out. Taking into account that the clus- 
ter volume is several cubic megaparsecs, the magnetic 
field provided directly by galactic winds and AGN lobes 
are grossly insufficient and the observed field has to be 
amplified in situ by cluster turbulence. Typically, the 
outer scales of cluster magnetic fields estimated from ob- 
servations (jCarilli fc Tav lor 2002) are much smaller than 
turbulence outer scales, which is why cluster dynamo is 
normally the small-scale dynamo. Also, generating large- 
scale fields, i.e. larger than an outer scale of turbulence, 
such as those observed in spiral galaxies, requires many 
(typically hundreds) dynamical timescales of turbulence 

^ A lower bounds based on the argument that FERMI did 
not observe GeV photons from inverse Compton scattering of 
the CMB by the electron-positro n pair cascade ha s been sug- 
gested by INeronov & Vovy ifSoToh : ITavecchio et al.l l|20Tol 'l. but 
has been deba ted recentl y by the beam instability argument by 
IBroderick"eral.i 120121 ): iSchlickeiser et al.l l iSOl^ ). However, an- 
other lower bound of SB = 10'^^ G{n / lO'^ cm-^){T / 10^ K)'^ / 
based on b a sic pl asma processes has been proposed recently by 
ISchlickeiseJ l[20ll) . 



on the driving scale, which is impossible in the ICM, 
where dynamical timescale is comparable with the age of 
the Universe. 

Although Coulomb mean-free path (mfp) of a thermal 
particle in the ICM is very large, 10-lOOkpc, it is be- 
lieved that the actual mean free path is much smaller 
due to magnetic fields and turbulence. Indeed, the Lar- 
mor radius of a thermal paricle in a 1/iG field is 10~^pc 
and even if we assume that particles stream freely along 
tangled magnetic field, the mfp will be greatly reduced. 
Further reduction of mfp is very likely due to scatter- 
ing by collective effects (iSchekochihin fc Cowley [20061 ; 
iLazarian fc Beresnvad 120061; iSchekochihin et al.l boOSD . 



with estimates of the mfp between 10~^ and 10~^pc. The 
cluster environment, therefore, is unlikely to be kinetic- 
viscous, as was commonly suggested before based on 
Coulomb mfp, but rather a high-Reynolds number {Re) 
turbulent environment, with self-similar scalings for ve- 
locity covering many orders of magnitude. While the 
kinematic dynamo, which ignored the backreaction of the 
magnetic field, has been studied extensively, due to the 
relative simplicity of the approach, the nonlinear small- 
scale dynamo received less attention. In high-i?e environ- 
ments kinematic dynamo saturates very quickly, giving 
way to the nonlinear regime. 

Nonlinear small-scale dynamo has been studied exten- 
sively only recently. In partucular, numerical simula- 
tions have established that the saturated state of such 
dyna mo is relatively unaff ected by Re, as long as Re is 
large (jHaugen et al.l[200^ . It was suggested that the 
small-scale dynamo in large Re flows could be univer- 
sal (Schliitcr fc Biermann 1950; Schekochihin & Cowle^ 
|2007[ ) and this was supported by further numerical 
and a nalytical studies (jCho et al.ll2009l ; iBeresnvak et al.l 
|2009| ). in particular the universality of the homogeneous 
small-scal e dynamo based on turbulence locality has been 
argued in iBeresnvakI (|2012l ). Although real clusters are 
not homogeneous, small-scale turbulence could be con- 
sidered approximately homogeneous on scales which are 
much smaller that cluster size. In this case the growth 
rate of magnetic energy is proportional to the local tur- 
bulence dissipatio n rate with a coefficient of around 0.05 
(|Beresnvak|[20Tl ). 

Full self-consistent MHD simulations of galaxy cluster 
formation, such as those presented here, are, in princi- 
ple, able to model the small-scale dynamo in clusters. 
However, due to the much lower effective Re, such sim- 
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ulations tend to prolong the initial exponential growth 
period. A balance have to be found between realistic 
magnetic injection mechanism, su ch as AGN activity in 
iXu et all (f2010) or galaxy winds in lDonnert et all ()2009l ) 
and the aforementioned artificial delay of gro wth. Since 
the ni agnetic fields obtained in simulations by iXu et aH 
((Ml) that we use in this paper are consistent with obser- 
vations, we will assume that they adequately reproduce 
the small-scale dynamo action in clusters. 

3. CLUSTER SIMULATIONS 

The ga laxy c luster studied in this paper is the simula- 
tion A in IXu et al. (2010). This simulation is performed 
using the cosmological MHD code with adaptive mes h 
refinement (AMR) ENZO-hMHD (|Collins et all IMoh . 
The simulation here uses an adiabatic equation of state, 
with the ratio of specific heat being 5/3, and does not in- 
clude heating and cooling physics or chemical reactions, 
which are not important in this research. 

The initial conditions o f the simulation are gen erated 
at redshift z = 30 from an iEisenstein fc Hul ()1999f ) power 
spectrum of density fluctuation in a ACDM universe with 
parameters h = 0.73, Vl^ = 0.27, = 0.044, = 0.73, 
and (78 — 0.77. These parameters are close to the values 
from WMAP3 observations (Spergel et al. 20Q2). The 
simulation volume is (256 h~^Mpc)^ , and it uses a 128^ 
root grid and 2 level nested static grids in the Lagrangian 
region where the cluster forms. This gives an effective 
root grid resolution of 512^ cells (~ 0.69 Mpc) and dark 
matter particle mass resolution of 1.07 x IO^'^Mq. Dur- 
ing the course of the simulation, 8 levels of refinements 
are allowed beyond the root grid, for a maximum spatial 
resolution of 7.8125 h^^ kpc. The AMR is applied only 
in a region of (~ 43 Mpc)'^ where the galaxy cluster forms 
near the center of the simulation domain. The AMR cri- 
teria in this simulation are follows. During the cluster 
formation but before the magnetic fields are injected, in 
addition to the density refinement, the refinement is con- 
trolled by baryon and dark matter density. After mag- 
netic field injections, all the regions where magnetic field 
strengths are higher than 5 x 10"® G are refined to the 
highest level. 

The ma gnetic field initializati on is using the same 
method in IXu et all (120081, 120091) as the original model 
proposed bv iLi et a l. (2006) assuming that the magnetic 
fields are from the outburst of an AGN. We have the mag- 
netic fields injected at redshift z = 3 in the most massive 
halo of 1.5 X 10^'^ Mq. We assume that the magnetic 
fields are from ^ 10^ M© supermassive black hole with 
about 1% of outburst energy in magnetic form. There is 
~ 1.9 X iqSO erg magnetic e nergy put into the ICM. Pre- 
vious study ()Xu et al.ll2010f) has shown that the injection 
redshifts and magnetic energy are not important to the 
distributions of the ICM magnetic fields at low redshifts. 

The simulated cluster is a massive cluster with its ba- 
sic properties at redshift z = as follows: Rviriai = 2.16 
Mpc, M^iHai (total) = 1.25 X lO^^ Mq, M„„ai{gas) = 
1.86 X lO^" Mq, and T„„ai = 7.65 keV. This cluster is 
already relaxed at the current epoch in the sense of X-ray 
and density distribution, but the turbulence is still ex- 
cited by recent minor mergers. The final total magnetic 
energy is 1.43 x 10^^ erg. The details about the cluster 
formation and magnetic field evolution are presented in 
IXu et al.l |2010, ). We visualized three quantities from a 




Fig. 2. — Second order structure functions, characterizing kinetic 
(solid) and magnetic (dashed) energy density in clusters. Here 
Pa, pb,Va,vi,, Ba, Bi, are the quantities taken at points a and b, 
separated by distance /. The middle point between a and b lies 
within a shell, whose number is indicated above the data. The 
outer radius of nth shell is n • 300 kpc. We use structure function 
method to calculate these quantities in shells around cluster center, 
since all quantities depend strongly on the distance to the center. 

simulation slice in Fig. [T] 

4. COMPRESSIBLE MHD TURBULENCE IN THE ICM 

These simulations reveal that the dominant mecha- 
nism of the excitation of the cluster-wide turbulence in 
the galaxy cluster medium is cluster mergers (see, e.g . , 
iVazza et al.|[20Tlt IXu et aLllMot iDonnert et aI|[20Tll) . 
Since the typical infall velocity and the typical thermal 
sound speed are related to virial velocity, they are of the 
same order, so the infall is typically trans-sonic, generat- 
ing a sizable amount of compressive perturbations and, 
possibly, weak shocks. Another mechanism could be due 
to bringing hot plasma inside cooler core environment, 
where the hot gas will become buyoant and produce con- 
vective turbulence and mixing in the center. 

Speaking of turbulence in the ICM, a simpli- 
fied approach is often taken when a turbulence 
is char acterized by the loc al cha racteristic ga s ve- 
locitv (llapichino fc Niemeveii 120081: IVazza et al.l 120091: 
IDonnert et al.l 1201 ID . This will not be sufficient, how- 
ever, if one wants to address the issue of second order ac- 
celeration. Indeed, the quasi-incompressible component 
of MHD turbulence, consisting of slow and Alfven modes 
tend to become progressively more anisotropic on smaller 
scales (jGoldreich fc Sridhad 119951 ) and, in a test parti- 
cle limit, effectively decouples from cos i nic ra ys (CRs) 
(IChandranI I2OOOI: lYan fc LazarTm] I200I. I200I . There- 
fore, we want to know the fraction of the fast mode, 
produced by cluster turbulence, which up till now has 
not been reliably estimated. In this paper we will es- 
timate the ratio of turbulent energy in each mode from 
simulation data and our subgrid model of turbulence will 
assume independ ent energy cascades for each mode as in 
ICho et al.l (I200IFI . 

The question of mode coupling is outside the scope of this 
paper, however a few comments could be made. In particular, 
the absence of energy exchange between slow and Alfenic mode 
on small scales in a weakly compressible case can be argued on 
rigorous basis, even if the interaction between modes is strong, i.e. 
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Fig. 3. — Turbulent dissipation rate in clusters, calculated with 
third order structure functions in three shells. Combined with 
the previous data, this allows us to estimate dynamic lifetime of 
turbulence in cluster, which is 3 ■ 10^ years in shells 2-3 and 8 ■ 10® 
years in shell 5. 

Estimating mode ratios is a fairly challenging task, be- 
cause the magnetic field in clusters is tangled and has 
no global mean field component, so global Fourier trans- 
forms of turbulent fields over the whole cluster could not 
be decomposed into modes. Also, doing global Fourier 
transforms in a datacube containing a whole cluster is 
rather meaningless, since all quantities strongly depend 
on the distance to the cluster center. So, the power 
spectra obtained from such transforms will be severely 
contaminated by large-scale gradients. In this paper we 
propose and use a local decomposition method based on 
structure functions (SFs). This is the first time such a 
method is used with the ICM MHD turbulence. This 
method also allows us to reduce the error associated 
with the fact that the cluster is not uniform and any 
measurement will include large-scale gradients associated 
with the cluster shape. In particular all global velocity 
measurements are contaminated by global accretion flow. 
The attempts to subtract this flow based on careful mod- 
elling of the averaged flow has not been either simple or 
reliable. 

The second order structure function is directly related 
to the power spectrum in the homogeneous case, e.g. the 
second order structure function scaling r™ corresponds 
to the power spectrum scaling of fc"™"^ (0 < m < 2). 
In the inhomogeneous case, the structure fuction will act 
as a good proxy for determining the local power spec- 
trum of turbulence. The main idea of our method is to 
use volume averaged structure function in shells centered 
on the cluster center. We are able to do so because the 
structure function is a local measurement. We used five 
shells with radii increasing linearly by 300 kpc for each 
shell. Fig. [5] shows second-order structure functions cor- 
responding to the kinetic and magnetic energy densities 
in the cluster. Note that the power spectra in different 
shells differ by orders of magnitute. This further reiter- 
ates the need to use local measurement, rather than the 
global Fourier transform. 

Fig. [3] shows turbulence energy flux, calculated locally 

perturbative theor y is im possible, see, e.g., IGoldreich &: Sridhari 
l|r995 ): Bcrcsnval3 l|20Tir i. 



by a well-kno wn method of third-order structure func- 
tions (see, e.g.. lPolitano fc Pouaueti[T998D . We can apply 
this method because cluster turbulence is only weakly 
compressive. Our energy fluxes correspond to the tur- 
bulence decay times of several billion years. This is a 
reasonable number for a quiet cluster. Note that if we 
assume that the cluster center is as hot as in our adia- 
batic simulations, the bremsstrahlung cooling in the cen- 
ter, ~ lO^^^erg cm~^s~^ will be a factor of ~ 30 larger 
than the turbulent heating rate in shell one. Turbulent 
heating is unable to compete with cooling. This is ex- 
pected, however, since it is turbulent mixing which is 
primarily responsible for heating the center. Indeed, if r 
is the turbulence mixing time, then turbulent dissipation 
is pii^/r, but the heating from mixing can be estimated 
as nT{l/lT)/T, where It is a temperature gradient scale 
and / is the turbulence outer scale. We expect turbu- 
lence to come from mergers and have outer scales as large 
as the cluster size, i.e. I /It ~ 1- On the other hand, 
nT ^ pv^ as clusters are subsonic. If Mg is an RMS 
sonic Mach number of cluster turbulence, the turbulent 
mixing provides appoximately 1 / times more heating 
of the cluster core than direct turbulent dissipation. 

5. MODE DECOMPOSITION IN MHD TURBULENCE 

MHD mode decomposition in turbulence, based 
on Fourier transforms ha s been use d before 

W [Cho & Lazariad (I2OOI : ICho et al.l (l2003h : 

iCho fc LazarianI (|20Q3| ). In these papers the global 
Fourier transform provided a global wavevector k for 
each Fourier mode, therefore necessitating the use of a 
global mean field. Indeed the global mean field was rela- 
tively strong in all aforementioned simulations, resulting 
in trans-Alfvenic turbulence. This approach works in 
the interstellar medium (ISM) of spiral galaxies where 
the turbulent component of the field is comparable 
with the mean- field component (see, e.g., iSchlickeiseii 
I2OO2D . In galaxy clusters such an approach is impossible 
due to the lack of a global mean field. The wavelet 
tec hnique has been pr o posed to deal with this difficulty 
in iKowal fc LazarianI ([2010'). The wavelet technique 
is computationally expensive, however, and requires 
datacubes with reduced resolution. Also, it has been 
only applied to homogeneous turbulence. The case of 
cluster turbulence is especially difficult, since not only 
cluster turbulence has no mean magnetic field, but also 
it is strongly inhomogeneous with density changing 
several orders of magnitude from the center to the 
outskirts of the cluster. In this paper we use a hybrid 
approach based on both Fourier transforms and struc- 
ture functions. As was described in previous sections, 
since structure functions are local measurements we 
can calculate them in shells around the cluster center, 
therefore mitigating effects associated with large scale 
gradient along the radius. 

As cluster turbulence exists in a very hot gas, ~ lO^K, 
the gas pressure is typically higher than dynamic pres- 
sure in most of the volume, resulting in subsonic turbu- 
lence with Ms <^ 1. Out of three MHD modes, Alfven 
mode is precisely incompressible, due to k being per- 
pendicular to both B and (5v. Slow and fast modes are 
both compressible, but the slow mode is almost incom- 
pressible in subsonic, high-/3 case. We can, therefore, 
split modes in two groups, one of them is almost incom- 
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Fig. 4. — Second order structure functions of Alfvenic (left), slow (middle) and fast (right) mode, calculated with respect to the local 
mean magnetic field, x axis is along magnetic field, while y axis is perpendicular to the field. Unit of structure function is lO^^cm-^/s^, unit 
of length on x and y axis is the box size, which is 3 Mpc. Alfvenic and slow mode eddies show small anisotropy along B, the anisotropy is 
small due to very short range of MHD scales, as the MHD scale is close to dissipation scale. Fast mode shows slight inverse anisotropy. 



pressible Alfven and slow mode and the other is fuUy 
compressible fast mode. In high-/3 case fast mode speed 
is close to sonic speed Cs and the dispersion relation is 
almost isotropic. 

Figure [4] shows the structure of the modes using con- 
tours of the second-order structure function plotted with 
X axis corresponding to the direction along magnetic 
field and y axis to perpendicular direction. The con- 
tors approximate the averaged shape of "turbulent ed- 
dies". We see that Alfvenic and slow modes show fairly 
small anisotropy along B, which is due to the fact that we 
have very short inertial range for MHD turbulence (below 
MHD scale). Fast mode shows slight inverse anisotropy. 
It was often assumed earlier that the fast mode turbu- 
lence is intrinsically isotropic. It is not true, however, 
since the weak cascade of the fast mode consists of inde- 
pendent cascades along rays in k-space, so it will preserve 
any anisotropy which is originally present. The scatter- 
ing by other modes is expected to be small on small scales 
and the cascade rate for each individual ray is a strong 
function of the angle between k and B. 

Figure [5] shows the second order structure function for 
velocity component corresponding to the fast mode. It 
turns out that the fraction e/ of kinetic energy residing 
in the fast mode is considerable. We studied this fraction 
by taking the ratio of the fast mode component second 
order structure function to the total velocity second order 
structure function and found that this ratio was around 
ef ^ 0.25 for all sh ells. This contrasts with the analyti- 
cal estimate used in iBrunetti fc LazarianI ()2011al ) , which 
gives the small fraction ~ M'^/Ma for subsonic clusters. 
The analytical estimate in the aforementioned paper was 
based on periodic box simulations with solenoidal driving 
that typically produces very small amount of fast mode 
(jCho fc LazarianI |2003( ) . In clusters, however, most of 
the turbulent energy is supplied due to mergers which 
have compressible and trans-sonic velocity fields, since 
the infall velocity is of the order of sonic speed which is 
virial on the outskirts of the clust er. This explains the 
differen ce between the estimate in IBrunetti &: LazarianI 
(j2011aD and our calculation of e/ from simulated cluster. 

Since the typical Mach numbers in the cluster center 
are modest Ms ^ 0.3, we will assume that the fast mode 
supplies a small fraction of the total dissipation rate, 



0.02eo- On the outskirts of the cluster this 
fraction is higher, due to Mg ^ 1, but the estimates of 
fast mode amplitudes on the outskirts are only tentative, 
due to limited numerical resolution. 

For the fast mode scaling we will use the so-c alled weak 
turbu lence model, which was suggested in (|Cho et al.l 
120031) . The power spectrum of fast mode fiuctuations, 
assuming their isotropy, which is roughly consistent with 
Fig. |5l can be expressed phenomenologically as 



(1) 



where Ck f is the Kolmogorov constant for weak turbu- 
lence. Note that the constant in this form is implicitly 
averaged over angle, since the weak cascade happens in- 
dependently along each ray in k-space. For the purpose 
of this paper, however, we will only need the assumption 
that the spectrum is approximately isotropic. 

6. IMPLICATIONS FOR CR ACCELERATION 

In this section we will briefiy outline the consequences 
of the significant amount of fast mode found in our clus- 
ter simulations. We will use simplified expressions for 
second-order acceleration, dropping out numerical fac- 
tors of order unity. More detailed calculation of acceler- 
ation, although based on the above mentioned analytical 
estimate of the amplitude of the fast mode, can be found 
in Brunctti & Lazarian (2011a, b). We also expect addi- 
tional factors from the fact that the fast mode, excited 
by cluster turbulence, is not exactly isotropic. 

Second-order acceleration by turbulence has been sug- 
gested as a process that provides additional energy to 
secondary electrons and allows them to produce an ob- 
servable amount of radio emission. Energetically, it is 
hard to support large amounts of fast electrons just by 
protons alone, and even if it was possible such sources 
will have much brig hter cores than observed (see, e.g., 
iDonnert et al.|[2"011[ ) . The attractiveness of second-order 
acceleration models is that it will produce relatively fea- 
tureless radio halos, similar to the observed halos, since 
turbulence is volume-filling. If, on the other hand, most 
of the acceleration happens in the accretion shocks, they 
will produce a different observed morphology, due to 
relatively short lifetimes of synchrotron electrons. In- 
deed, so-called radio relics, which morphologically re- 
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Fig. 5. — Second order structure function of the fast mode com- 
ponent of velocity, where the vector between points a and b was 
chosen to be perpendicular to the cluster center to exclude contri- 
bution from the accretion flow. Note that due to the poor resolu- 
tion at the outskirts of the cluster (shells 4-5) the scaling is very 
steep, i.e. turbulent scaling is not observed. This allows us to es- 
timate the amplitude of the fast mode perturbations only in the 
inside shells (1-3). 



semble l arge-scale shocks, h as been observed in some 
clusters (jEnsslin et al.|[l998[) . Although we can not ex- 
clude shock acceleration in ordinary radio haloes, due 
to either face-o n shock geonietry or mu ltiple small-scale 
shocks (see, e.g. Ivan Weeren et al.ll201Clt ). in this paper we 
will investigate primarily the volume-filling second-order 
acceleration, which is the most natural explanation of 
megaparsec-scale radio haloes. 

Second-order acceleration is a diffusion of particles in 
the momentum space due to collision with magnetic 
irregu l arities . The original mechanism, proposed by 
iFermil (|1949l ). assumed collisions with large-scale clouds, 
while turbulence models assume resonant interaction 
with MHD modes. The flux of particles through momen- 
tum space due to momentum diffusion is defined as = 
—'iiTp^D.ppdf/dp, where we assumed isotropy and have 
used the momentum diffusion coefficient Dpp already in- 
tegrated ow£lh£_2itch_arigle. This Dpp is equivalent to 
A2 from iSchlickeised (|2002 ) . This term provides acceler- 
ation, as long as df /dp < 0, which is normally satisfied. 
The estimate of momentum diffusion due to the interac- 
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tion with fast mode as Dp^ ^ ^ ... 

where Lm is the MHD scale (see, e.g.. lSchlickeisedl2002l : 
lYan fc LazarianI [20041 ). We dropped out all factors of 
order unity, although they can also be absorbed in the 
parameter Lm- If we assume that this process is domi- 
nant over space diffusion and losses, then the stationary 
solution will correspond to = const which suggests 
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corresponding to the energy distribu- 



tion Ep — 4,11 p^ f P~^l'^^ a very flat spectral index, 
compared with observed spectral indices which are be- 
tween —2 and —4. This flatness is an inevitable feature 
of second-order acceleration, because particles try to dif- 
fuse to higher p where there is a larger phase volume. 
If escape due to diffusion is taken into account, the re- 
sult is virtually unchanged as long as particles are well- 
trapped, i.e. acceleration timescales are smaller than the 
diffusive escape timescales. If, on the other hand, escape 
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Fig. 6. — Spectra of electrons and protons, reaccelerated by fast 
mode turbulence in the galaxy cluster. The proton spectrum is 
limited from above by absolute energetic constraint of Eq. l[2)l. 
The e/p ratio initially corresponds to the injection e/p ratio till 
the electron synchrotron suppression at ps- This ratio will increase 
later due to secondary production. For higher energies e/p 
ratio is smaller, as electron spectrum is in equilibrium between 
production directly from protons and synchrotron losses. 

is more efficient, the second-order acceleration will only 
slightly modify th e injected particle distribution (see, e.g. 
iSchlickeiseii 120021 ) . The radio emission spectra of halos 
are rather steep, which suggests steep particle distribu- 
tions and it would seem that second-order acceleration 
will either be inefficient or produce particle distributions 
which are clearly incompatible with the observed ones. 
However, in practice a variety of other effects have to 
taken into account. Synchrotron cooling is one of these 
effects, it could be expressed as = —Anp^pf, where 
p ^ p^. This will result in a cutoff of accelerated distri- 
bution, in a form Ep exp(— (p/ps)"^/^), where p^ 
is a cufoff momentum. This cutoff will result in a steeper 
emitted radio spectra and these spectra can be mistak- 
enly identified as produced by steep Ep distributions. 

We argue, however, that another process of fundamen- 
tal importance is operating in second-order acceleration, 
which is the back-reaction of particles to turbulence. In- 
deed, there is an energetic constrain that requires that 
CRs should not extract more energy than that is avail- 
able in the turbulent cascade. If this condition is ignored 
and turbulent spectra are postulated to be power-law, 
an excessive amount of particle heating will be obtained. 
The turbulent cascade can produce classic self-similar so- 
lutions, such as SB^ ^ k^^^^ for the fast mode, only 
when the back-reaction to turbulence is ignored and the 
energy flux is due to fluid nonlinearity alone and is con- 
stant through all scales. In practice, this condition will 
be quickly broken by the feedback of CRs. 

An absolute upper limit on second-order acceleration 
could be obtained by assuming that the energy extracted 
on each scale should not exceed large scale turbulent dis- 
sipation rate e. As the flux of energy in accelerated parti- 
cles can be expressed as Fe — pcF — —Aircp^Dppdf /dp, 
this upper limit will result in an upper limit on Ep < 
p~^'^, which is steeper than constant particle flux solu- 
tion, suggesting that unless second-order acceleration is 
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limited by losses or escape, it will eventually be limited 
at high energies by back-reaction to turbulence. Assum- 
ing that protons take most of the turbulent energy the 
Ep will be limited from above by 

Ep < ev-/LMP^^'^p-^'\ (2) 

where Lm is the MHD scale of fast mode turbulence, 
Po « 12eBLM / c. The transition to steeper spectra pos- 
sibly generated by other mechanisms happens at a mo- 
mentum pd, which will be determined by the total turbu- 
lent energetic history of the cluster and is approximately 
related to the temperature on the outskirts of the cluster. 
Assuming turbulence pumping 30% of the time with the 
dissipation rate of e during the age of the Universe i^, we 
obtain the expression pd — pii{Q.lbtAv\/cLMY- Numer- 
ically it could be estimated assuming va ^ 10^ cm/s and 
Lm ^ 60 kpc as 2 • 10^% - 2 • lO^^ eV/c, so the 
maximum proton gyroradius is still much smaller than 
MHD scale Lm- This is due to second-order acceleration 
being a rather inefficient acceleration mechanism. 

As for the electrons, their cutoff energy ps 
will be determined by the balance between syn- 
chrotron cooling and momentum diffusion as Ps — 
mec{lQTT'^eBm^c^ / a'^LMBlf/^ivAl cf'^ , where = 
(B2 + _b2^)1/2 includes equivalent B/c 3(H-z)2^G to 
account for inverse Compton losses on CMB and syn- 
chrotron losses. Using parameters from above, Ps 
4 • lO^TTieC. The sketch of the electron and proton spec- 
tra are presented in Fig. [6l The sharp cutoff in electron 
spectra, associated with transition from the spectrum of 
accelerated electrons to the spectrum of secondary elec- 
trons could explain the relatively steep spectra of radio 
emission from radio halos. 



7. SUMMARY 

The nonthermal radiation from hot conductive intra- 
cluster medium (ICM) is explained by the presence of 
accelerated particles. It is plausible that the cause of ac- 
celeration, in the case of radio halos, is not large-scale 
accretion shocks, but rather volumetric turbulence that 
is produced by constant mergers and buoyancy of the 
cluster medium. This turbulence is currently being in- 
vestigated by cosmological MHD simulations as being a 
source of both magnetic fields and accelerated particles. 
The key to in-situ acceleration of electrons and protons 
is the MHD fast mode which effectively scatter parti- 
cles and diffusively accelerate them. We propose a novel 
method to extract fast mode from MHD simulation and 
estimate the efficiency of acceleration and scattering. An 
important limit, imposed by the energy flux constraint 
of the fast mode turbulence leads to a spectrum which 
is steeper than one would naively assume from diffusive 
acceleration. Moreover, this spectrum has an exponen- 
tial cutoff at relatively low energies for electrons which 
could naturally explain very soft radio spectra typically 
observed for diffuse radio halos. A more detailed calcu- 
lation of the particle distribution in space and simulated 
radio maps will be presented in a future publication. 
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